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We study the contribution of two-body meson-exchange processes to axial 
charge transitions for nuclei in the lead, tin and oxygen regions. We conduct 
calculations in the Dirac-Hartree (the Walecka model) and the relativistic 
Hartree (where the full one- nucleon- loop effects are included) approximations. 

X' 

■ We present results indicating that one- and two-body processes enhance the 



matrix elements of the axial-charge operator by some (100ib20) % in all three 
regions studied. This agrees well with the fit of eighteen first-forbidden /3- 
decay transitions conducted by Warburton in the lead region. We also discuss 
some sensitivities present in the calculation. 

I. INTRODUCTION 

The unexpected enhancement found in processes driven by the axial charge, namely 
the cases of first-forbidden beta decays in heavy nuclei 0] and the production of pions near 
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threshold in pp — > ppir^ reactions [0], has given a new impetus to the study of meson-exchange 
currents (MEC) in axial charge transitions. The findings are really remarkable. In the first 
process for example — already known to be strongly renormalized by the presence of MEC — 
the results found by Warburton exceeded the anticipated theoretical estimate by a factor 
of two. In the second — a new result in a region not available to experiment before — the data 
overestimate the theoretical predictions (only the one-body contribution) by approximately 
a factor of five. Despite this sizable difference between theory and experiment the first 
result is acknowledged as the most dramatic since the process of first-forbidden /3-decay was 
thought to be (at least to a large extent) under theoretical control 0. 

The subject of axial charge transitions becomes more relevant from the viewpoint of 
relativistic theories of the nucleus. Most of the explanations advanced so far to account for 
these enhancements (refs. [^,0 in the first case and in the second) rely on the presence 
of scalar heavy-meson exchange between two nucleons. Diagrams containing a nucleon- 
antinucleon pair originating (or decaying) into a scalar meson of rris ^ 500 MeV using 
couplings from either the Bonn or Paris potentials seem to take care of the discrepancy 
between the theoretical prediction and the experimental result. These large Lorentz scalar 
potentials have long been a distinctive feature of relativistic calculations whose explicit 
need, however, had not shown up in the non-relativistic phenomenology of nuclear processes 
before. The calculations alluded to above (non relativistic) incorporate the scalar degree of 
freedom formally via meson-exchange currents. 

The connection between non-relativistic pair-exchange currents and relativistic mean- 
field theories for isovector currents was studied by Delorme and Towner in their paper 
of 1987 (in the case of isoscalar currents there were previous calculations by Arima et al. Q 
and Blunden ||^). Pair graphs are shown in Fig.|l| (diagrams (a) and (b)). The relativistic 
Hartree contribution is obtained by closing the line on nucleon 2 and is shown in diagrams 
(c) and (d) (diagram (d) is usually recognized as the "backfiow" correction). The difference 
between both approaches is that whereas in the relativistic case the propagation of an 
antinucleon is included in the one-body Green's function, in the non-relativistic case it can 
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only be incorporated as an exchange current. Delorme and Towner found in their calculations 
large enhancements due to the scalar-exchange terms. Their conclusion at the time — for 
axial-charge transitions — was that such enhancement (and along with it the relativistic 
mean-field theory) was hardly consistent with experiment. Quite oppositely, they noticed 
that the conventional theory of meson exchange with short range cut-offs gave reasonable 
agreements with the data (ignoring, of course, the contribution of pair diagrams). The more 
recent work of Warburton on the one hand and the results of Meyer et al. on the other, 
prompted a reanalysis of the existing calculations under the light of the new data. At the 
present time it is generally held that scalar heavy-meson exchanges can account for the 
discrepancies between the predictions and the actual data [ p!0[] . 

This situation opened up the ground for further exploration of the predictions of rel- 
ativistic models of nuclear structure to axial-charge processes. Since the physical content 
of the two descriptions (non- relativistic vs. relativistic) is different from the outset, the 
individual contribution of one- and two-body diagrams is not anticipated to be similar. We 
already mentioned that the one-body relativistic calculation includes the non-relativistic 
two-body pair graphs. Conversely, the two-body contribution in the relativistic calculation 
sums diagrams (like the one shown in Fig.|^) that are not included in the traditional meson- 
exchange calculations. Though, a priori, contributions of this type are small a comparision 
between both approaches will only be meaningful once full calculations of one and two-body 
processes are conducted within each framework. 

With this in mind we study in this paper the contribution of two-body processes within 
the relativistic approach. We concentrate exclusively on first-forbidden beta decays in nuclei 
(oxygen, tin and lead) where the relativistic mean-field can be used reliably. Regardless of 
the model the symmetries must be upheld. In particular chiral symmetry, namely the use 
of soft-pion theorems and current algebra, imposes a strong renormalization of the matrix 
elements due to pion-exchange processes (particularly through the vrp-exchange diagram.) 
This mechanism needs to be incorporated explicitly into relativistic calculations since, at 
the mean-field level, the pion is absent due to parity arguments. It is clear that a relativistic 
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description of axial-charge transitions in nuclei requires to go beyond the mean-field. A first 
approximation to the problem calls for a mean-field calculation (impulse approximation) 
to describe the one-body processes and additionally a MEC calculation taking care of the 
"backfiow" and pion-exchange processes. We notice that the results of this programme are 
not expected to be free from the usual ambiguities of pion couplings in relativistic theories 
(for a review of the current status see ref. |lTl[| ). 

Recently, within the framework of a relativistic model (quantum hadrodynamics II) 
we analyzed axial transitions |]12| in heavy nuclei using the scheme outlined above. We 



obtained first the one-body contribution given by graphs (c) in Fig. |l| next, the RPA- 
type vertex corrections (graphs (d) in the same figure) ||13[ and, finally we investigated the 
contribution of the pion-exchange currents 0] (with no exchange diagrams.) Inspection of 
the results arrived at in this study showed that: First, despite some ambiguities related to 
the strength of the pN coupling constant, the sum of one- and two-body diagrams added up 
to an enhancement which was comparable to that obtained in non-relativistic calculations. 
Second, the theoretical calculations could only be brought into agreement with the data 
when pseudovector vrN coupling was used (otherwise the results were anomalously large). 
Third, two-body pion-exchange processes were indispensable to understand the axial-charge 
data. These conclusions were based on calculations conducted on two particular transitions, 
one in Tl and the other in Pb, and were intended to give an idea of what should be expected 
for other transitions in this mass region. In this paper we extend and deepen the analysis 
of the previous works. On the one hand, we give a complete description of the two-body 
pion-exchange contribution to the axial charge from the point of view of relativistic hadron 



theories. On the other, we extend the results of |jT^ and present results of MEC contribution 
to axial charge matrix elements for tin and oxygen. The idea is to move to lighter nuclei to 
study mass dependence, and the infiuence of the scalar field in this contribution. 

Besides this introduction the paper contains two other sections. Section divided into 
two subsections, deals with the calculation of the matrix elements of the axial charge. In 
the first subsection we give some details about the hadronic model we employ and the way 



the axial two-body current is derived. In the second we describe at some length the scheme 
to carry out the actual calculation of the matrix elements of the current between the nuclear 
states. The last section contains a discussion of the results and of the potential ambiguities 
involved in the calculation. Finally we present some conclusions. 



II. MESON-EXCHANGE-CURRENT MATRIX ELEMENTS OF THE AXIAL 

CHARGE 

A. Hadron model and currents 

At tree level the two-body isovector axial vector current is given by the sum of the 
four diagrams shown in Fig. |^. In the figure we have used the symbol e to denote the field 
operator mediating the weak current. Diagrams a) and h) turn out to be the most important 
since they involve pion exchanges which generate the longest-range exchange currents. The 
contribution of the other two diagrams (involving the A) is known to be small [T3|. 



The two-body current from diagrams a) and h) can be calculated using standard Feynman 
rules. As stated meson and baryon propagators derived from some baryonic Lagrangian are 
needed. Throughout the paper we use the quantum hadrodynamical model 11 (QHD-II) of 
Serot and Walecka . It includes the nucleon fields, two isoscalar meson fields -the scalar 



a and the vector uj- and two charged isovector meson fields -the pseudoscalar vr and the 
vector p. The full Lagrangian expression is, 

+gpid''n X 7f) • p; + ^gliTT x p^) ■ (vf x p^) + ^{d^ad^a - m'y) (2.1) 

+ ^((9^7? ■ 9^7? - m^vr ■ 7?) + ^g^^m^an ■ vr 
In Eq. (|2.1| ) the neutral and charged vector meson field strengths are defined by, 

F'"' = d^uj" - d^uf" (2.2) 



and 

B^,u = df,pu - duP^, - Qpip^, X p^) (2.3) 

respectively. We have also made use of the covariant derivative defined as 

1 

D^ = d^ + i-gpf- p^. (2.4) 

There is a Higgs-meson sector associated with the p-meson field which has been suppressed 
by taking the scalar Higgs mass to be very large. 

The ttN coupling in the QHD-II Lagrangian is pseudoscalar. To account for the small 
pion scattering length in vrA^ interactions a crTrvr coupling term is added. As we will see 
below pseudoscalar coupling with the standard parameters gives very large predictions for 
the contribution of one- and two-body processes. As an alternative to the actual calculation 
we have considered substituting pseudovector for both pseudoscalar vrA^ and cttttt couplings. 
We emphasize that despite the different predictions both couplings show the same correct 
behaviour in the low-energy limit (in the two-body case) and both satisfy the constraint of 
partial conservation of the axial current (PCAC). 

The strong vertices ttNN, pNN and pvrvr as derived from the Lagrangian of QHD-II are 
shown in Fig. The weak vertices ne y Tcpe (see Fig. H) cannot obtain directly from the 
QHD-II Lagrangian since it does not contain the weak interaction. The structure of the ne 
vertex is obtained from symmetry arguments in pion decay where, 

r^f = {vac\ J^'^" I n\q) ) = (2.5) 



The value of the constant, = 0.93 GeV, is taken from experiment |T8|. For the npe vertex 
a similar argument leads to a vertex of the form, 

^,uabc ^ ^ ^a^f^^-^ ^6(^^^ | jA,,c | ^^^^^^^ ^ ^^^^ ^3.6) 

where now the structure and content of the T^'^ tensor cannot be determined from Lorentz 
covariance only. However we will see below that the requirement of PCAC on the two-body 
axial current suffices to define it. 



The propagators of the vr and p fields appearing in QHD-II are of the Klein-Gordon type. 
With them and the vertices just defined we can obtain expressions for the contributions to 
the two-body axial vector current density following Feynman prescriptions for diagrams a) 
and h) in Fig. ^. For the case of pseudovector vrA^ coupling they read, 

J^i"' = -iff In X T,r *(p;)tS-a*;?*(pi) 

for diagram a) and, 



AM k( — 

xg- (Q + hY -^^^^W2)1uHp2). (2.8) 
q — K2 — iiip 

for diagram b). 

The sum of ( p.7| ) and ( |2.8| ) gives the full two-body (2B) isovector axial-vector current. 



X i^T^' + 9ph^^ {q + hY } i^(P2) 7. ^(P2) (2.9) 

The quantity between brackets behaves as the effective irpe vertex with the tensor T^^ still 
an unknown. This last difficulty can be overcome after enforcing the requirement of PCAC 
on two nucleons. This is equivalent to demand that the four-divergence of the full two-body 
current be equal to, 

2 

q J^.M^ = , M%2B) (2.10) 



In ( |2.10|) , M^{2B) is the vr-absorption amplitude on two nucleons which obtains from the 



pixTi diagram (diagram h) in Fig. ^ by removing the external e^^^"" line. Thus from the figure, 
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X + T^—^ T^—-, ^iP2) lu^{P2). (2.11) 
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Substitution of ( |2.11| ) into the right hand side of ( p.lO| ) gives for the T^'^ tensor the expression, 



rpflU 



'fngp 



q2 _ ^2 



(2.12) 



thus fixing the strength at the pvre vertex. It is important to notice that in the soft-pion 
hmit {ki —>■ 0), one obtains T^^" = —fngpg^^- This is the same strength obtained by taking 
the soft-pion hmit in hard-pion models (see for example, Ivanov and Truhlik [pj|] .) 

Substituting T^^ into the expression for the axial current and using the relation fj^g^^ = 
qaM (Goldberger-Treiman), 



J. 



A,fi c 



2B 



.9A9f: 
I 
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1 



kf — ml k^ 



q2 _ ^2 



(2.13) 



which is our final expression for the full two-body isovector axial- vector current. Inspection 

of the quantity between curly brackets shows that J^'^'^ satisfies PC AC since, 

,2 



jA,ric I 

Qp J2B \ -9 



q2 _ ^2 



mz 



^ 



(2.14) 



Thus, it is clear that the two-body axial current derived from the Lagrangian ( p.l| ) satisfies 
PC AC and has the correct behaviour in the soft-pion limit. 



B. Multipole decomposition 

Having studied the structure of the relativistic two-body current we concentrate next 
on the nuclear structure. Individual transitions in nuclei are more conveniently studied 
by performing a multipole decomposition of the current matrix elements. For axial-charge 
mediated processes selection rules limit the contributing operators to the axial-Coulomb 
(or charge) operator, and the polar-vector longitudinal and transverse-electric operators. 
The last two affect only minimally the total matrix element and therefore, though they will 
be taken into account in the actual calculations and displayed in the results, for the purpose 
of this section we shall work exclusively with the first one. The axial-Coulomb operator is 
written in the form, 
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Ml 



t2B, 



X) 



(2.15) 



with 



Mf{qx)=Mqx)Yr{x) 



Ml 



(2.16) 



The reduced two-body matrix element of any multipole operator may be decomposed 
as a sum of products of a matrix element of the operator between two-particle states times 
a two-body density matrix element which accounts for the complexity of the many-body 
problem. Thus, for the axial Coulomb operator we shall write, 



{Jf\\tT\q)U)=Y. E <[a\h']j,\\6r\maMji>'^'}'''{[<^\h']^^^^^^ (2.17) 

a,b,ji a',b',j2 

where the sums run over all initial (final) two-particle states coupled to angular momentum 
ji (j2)- The two-body density matrix element is written in terms of single particle creation 
operators, cj^, in the form 



¥i'^\[a',b']n,[a,b]j,)) 



[J] 



< Ji 



c{, ® c{, 



32 



[Ca ® Cb]^^ 



\Jf > 



(2.18) 



with Ca = (—1 )•'""'"'"'' c_a. (We refer the reader to Sec. IV. E in Donnelly and Sick's review 
article |^ for the explicit expressions of the 2B-density matrix elements in the case of one 
particle or one hole outside a closed shell which are those of interest for this work.) 

Equation ( p. 17 ) requires the calculation of the reduced matrix element of the axial 
Coulomb operator between two-particle states, 



[a',b']Jl')= J2 {jamajbmb\jimi)\jama)\jbmb) 



(2.19) 



ma,mi, 



where 



{X\jama) = Kki^) 
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exp{-iEnjJ} 



(2.20) 



with 

X = ^ < Imlmsljfi > rrxT (2.21) 

and 

K = (j„ + (2.22) 

Since single particle states are obtained in co-ordinate space we transform the current ac- 
cordingly. 



1. Coordinate space two-body matrix elements 



Our starting point will be the total 2B current of Eg. (|2.9|) (where, for convenience, all 
but the Lorentz indeces have been dropped.) In the Heisenberg representation we write the 
current in co-ordinate space in the following form, 

1 f d^^d 1 



X J d'x, Alix - X,) dl [7/;(xi)7V^V(a:i)" 

X / d'x2 A%x - X2) h^?^"] \^Pix2h,T'i^ix2)] (2.23) 



where 



^„^_|-,-,.,.-0^_ ,,,,, 



and M'fp = +^ml. ^ — (n^y. For the pNN vertex we have worked in the t'Hooft-Feynman 
gauge. We have kept a soft q-dependence in the second term of the npe vertex in (|2.23|) that 
will be dropped in actual calculations. For future use we expand the meson propagators into 
multipoles using, 

p—m\x—x'\ oo 

. I ^, I = -mY^ ji{imrU)hf\imr' ^)Y*^Xn)Yi,^X^') (2.25) 
An \ X — X \ , 

I I l,mi 

where ji{z), and hf'\z) are spherical Bessel functions, and r'<(r'>) denotes the lesser 
(greater) of | x | and \ x' \. 
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To calculate the reduced matrix element < [a' ,b']j2\\C^j^\q)\\[a,b]ji > we substitute the 
expansion ( 2.25 ) into ( 2.23|) perform the time and angular integrations and contract the 
Heisenberg operator between the appropiate two-particle states. For convenience we divide 
the total calculation into two separate contributions stemming from each of the two terms in 
the irpe vertex. For the first term proportional to the time component of the metric tensor 
(labelled A) we arrive at an expression of the form, 

< [a',b']j,\\Cjiqma,b]j, >= ^g^gj V2JTl{-iy^' J drr'jjiqr) 



X 



X <^ {Ea' - Ea)[Ga'{n)Fa{n) - Fa\n)Ga{ri)] 



+ 



GAr.) \ir + y] Gain) + FAn) |; + ^ 1 F.{n) 



^2=0 



x[GAr2)Gt{r2) + FAr2)F,{r2)] 







3a' Jb' 32 
ja jb Jl 
k h J 



where M^'p" = Jml - {Ea^ - Ea'^'Y, and 



lAk\ k) = ^ V(2j + 1)(2/ + 1)(2L + 1) [l + (-1)'+' +^ 



-T 



(2.26) 



(2.27) 



V 2 ^ "2 / 

The symbol ^± in ( [^.26D denotes the matrix element of the charge-raising (or lowering) 
isospin current between two-particle states, 



^± = < ci'b'l — [n X T2]^. \ab>= ± < a'b'\{'riT2 - T^T2)\ab > 



± 



i^a'p^ap — ^a'n^an) ^b'P^b'^ " ^a' ^A(ip (^b' p ^bp — h'n^bn) 



(2.28) 



The second term of the axial-charge operator proportional to qQdtj/{q^ — m^) contracts 
with the matrix in the pNN vertex giving rise to a time-like and a space-like contribution. 
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The first one is simply a renormalization of the matrix element of the operator Cj just 
calculated and that we now redefine so that, 



< [a\b']j2\\Cf{q)\\[a,b]n >=< [a\b']j2\\Cj{q)\\[a,b]n > ( 1 - 



— ml 



the second (Cf ) results in the following expression, 



<K6']j2||Cf(g)||[a,6]ji >=-—g^gl 



u i-iy I 



dr 



X jdnrlM^'^f^ ^2h + l3i,{iM:^T\)hi,{iM^^Tl)Ii,{K', -K) 



X - E,)[GAn)F,{n) - F,,{n)Ga{n)] 



+ 



X 



«2=0 



oo 

X ^(2/ + 1) 



/=0 ^(Z + Z2 + 2)(/ + Z2 + l)(^ + ^2) 

X [G';,.(r2)F,(r2)//i,i,(A;;.,, -A;^) - F,,(r2)G'6(r2)//M2(-^6', h)\ 



X 



^/2ili^{kv,-h)Zj{iMM)< 



3a' 3b' 32 
3a 3b jl 
ll 12 J 



where the functions Hi^i^{ky, kb), I^i^ikb', kb) and Zj{l, li, I2) are: 
Hi,i^{kb', kb) = 5ii^\/2{ky - kb) + 5i^i, 



(2.29) 



(2.30) 



(2.31) 



W^b',h) 



40F 



2jb, + l)(2j, + l)(2/2 + 1) 1 + (-1)'"'+'" 



/ • / ■ \ 



1 n 1 

V 2 ^ ~2 y 



(2.32) 
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+ Jj{2J-l)jj_^{qr) 




(2.33) 



Finally, the sum of ( p^.29D and ( |2.3(J| ) gives the two-body matrix element of the axial- 
charge operator. 



III. RESULTS AND DISCUSSION 

The coupling constants used in the calculations are shown in Table |. Both sets of 
parameters, DH (Dirac-Hartree) and RHA (relativistic Hartree approximation), are from 
Ref. and are not tuned specifically to any particular mass region. The scalar masses 
were determined by fitting the charge radius of ^°Ca. The other masses, fixed to their 
experimental values are, M = 939 MeV, = 783 MeV, rup = 770 MeV and = 139 
MeV. In the calculations we used from the fitting to the bulk symmetry energy 04 = 35 
MeV. The intrinsic structure of the nucleons was considered via an additional dependence 
of the form (1 + gV^l)"^ with Ma = 1.2 GeV. 

Tables to summarize our results. In the tables each row displays the result of the 
calculation of the matrix element corresponding to the multipole operator shown on the left. 
They are the axial Coulomb operator (Eq. |2.15|) and the polar-vector longitudinal 



1 d^x 


-VMf{qx) 







(3.1) 



and transverse electric, 



-V X Mf (gr) 



J'^(x), 



(3.2) 



operators. In the above we used the expression 



(3.3) 
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with 

y!i = Y.c'^'ML{^y,- (3-4) 

the vector spherical harmonics. 

Each table is divided in three columns: The first one (NR) shows the results of a one- 
body non-relativistic calculation and has been included, mainly, for comparative purposes. 
The second (DH) shows the results of: i) a one-body (lA) calculation, ii) one- plus two-body 
(lA+MEC) calculation and iii) the resulting increment (5a), all three obtained within the 
relativistic Dirac-Hatree approximation where the vacuum structure is ignored. The third 
column has a structure similar to the previous one but the one-loop vacuum corrections that 



originate due to the presence of the strong meson fields are taken into account p2 . 

The increment 6a defined in the tables parametrizes the difference between the "effective" 
axial coupling constant, as determined from the analysis of the experimental transition rates, 
and the "free" coupling constant such that, 

=9^r{^+^A). (3.5) 

We recall that the "best fit" value reported by Warburton in his analysis of nuclei in the 
A=206-212 mass region was 

5a = 1.05 ± 0.05. 

The first thing to notice from the tables is that, as indicated above, the longitudinal 
and transverse electric contributions in the one-body case are substantially smaller than 
their Coulomb counterparts (more than a factor of ten in lead and tin; four times smaller in 
thallium). Based on this result we restrict the two-body calculation to the contribution of 
the axial charge where, on the other hand, we expect the larger contributions coming from 
pion exchange. 

Despite not being the main thrust of this paper we devote one paragraph to comment 
on the one-body vector contribution where we notice the following. In the lead transition 
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(table IH), the lA result for and is enhanced a barely 3% in DH (1% in RHA) with 
respect to the NRIA calculation, a result that is consistent with Warburton's finding for the 
polar vector enhancement, namely 5y ~ 0. This behaviour does not maintain, however, and 
we find that in the thallium decay to lead the enhancement is 65% (34% in RHA) and in tin 
is 23% (18% in RHA) over the NRIA. Thus, there is a discrepancy between this calculation 
and the fit to the data in the heavy-mass region |jl|| for the vector contribution. A priori 
this difference could be adscribed to the simple-minded model used here to describe the 
structure of the states involved; however is hard to see how a refinement in this direction 
could reduce drastically the last two values. 

Back to the analysis of the tables the next thing we notice is that, as expected, the one- 
body calculation is consistently larger in the Dirac-Hartree approach than in the relativistic 
Hartree approximation. This is due to the value of the effective mass M* which is smaller 
in RHA than in DH. In principle, there is no theoretical argument favoring the use of one to 
the other given that, in both cases, the model parameters are fitted independently to bulk 
nuclear-matter properties. For convenience, in the analysis that follows we shall limit our 
discussion to the RHA results. 

In the Introduction we argued that in the relativistic framework pair diagrams are already 
taken into account in lA; hence the analysis must proceed by studying the contribution of 
MEC along with the lA results. Thus, the u : 2gg/2 ^ tt : 1/19/2 transition in 209p]^ _^ 209gj 
shows that some 40% of the enhancement is present at the one-body level and a further 80% 
shows up at the two-body level bringing their sum up to a total enhancement of 6a = 1-2. 
The other transition we studied in a heavy nucleus, tt : (3si/2)^^ ^ '■ (3pi/2)"^ in ^'''^Tl — *■ 
207p^ gives a 30% enhancement due to one-body processes and 55% coming from the meson 
exchange contributions. 

If we move to lighter nuclei the situation does not change much as far as the overall incre- 
ment of the effective coupling constant is concerned. In the ^'^^Sn — i^^^s decay, dominated 
by the single particle transition v : 2/7/2 — > vr : l(?7/2, we obtained 42% enhancement from 
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lA and 33% from the MEC processes. In the oxygen region the decay ^^N(0~) — 0(0"'") 
was studied using a simple model of nuclear structure where the first excited state of ^^N, 
J'^ = 0" (with an excitation energy of 120 keV) was described by the mixing of two particle- 
hole configurations, 



I 0- ) = V1-A2 I [2si/2lpi/2]o- ) + A I [lc/3/2lp3/2]o- ) (3-6) 
and where the mixing amplitude A was taken to be the typical A ~ 0.08 ||23|. In this case 



we obtained 48% (lA) and 42% (MEC) for the largest component, and 18% (lA) and 10% 
(MEC) for the smallest. 

Overall the calculated value of Sa is around 20% of that determined empirically in the 
heavy-mass region. The ratio of one- to two-body processes, however, is not distributed 
uniformly in mass. As one moves towards the heavier nuclei the share between the one-body 
and two-body processes is shifted from the former to the latter. It is not straightforward to 
explain this trade off. Naively, it could be argued that the increase in the density implies 
a reduced effective mass and a corresponding enhancement in the nucleon propagators; one 
such change, however, should affect both the one-body and two-body processes in roughly 
the same manner. Furthermore, the argument of enhanced matrix elements due to off- 
diagonal operators (operators connecting upper to lower components of the spinors) it is 
not strictly true for /3-decay transitions where the initial and final states do not coincide. 
The value of the matrix element depends on the overlap of initial and final wave functions 
and this overlap diminishes as we go up in mass in the periodic table. A case in point are 
the different enhancements obtained for the vector operators, and ^\ which are very 
different in different transitions despite being both off-diagonal. Summarizing, the results 
of the tables indicate a somewhat stable contribution of the one-body processes of around 
40% whereas the two-body share seems to increase slightly with the mass (0:42%; Sn:33%; 
Tl:55%; Bi:80%). 

The calculation in oxygen with a simple nuclear model gives an enhancement close to 



100%. Before the work of Haxton and Johnson in 1990, this result would have been 
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considered anomalous. These authors found that the usually accepted 50% enhancement 
of the coupling constant fell short by some 40% of the value necessary to agree with their 
calculation of the /5-decay rate of the 120 keV 0~ state in nitrogen. This discrepancy plus 
the one found in heavy nuclei prompted a more careful examination of previous calculations. 



Presently, this discrepancy is not considered dramatic in any mass region [|IU| though there 
seems to be still some room for other effects. Our results seem to indicate a similar behaviour 
(so long as we stick to the RHA model and use pseudovector coupling). 

Regarding the type of coupling, Table ^ illustrates the difference between pseudoscalar 
and pseudovector contributions for a particular transition in lead. Use of the former leads, 
typically, to an increment of the calculated enhancement of approximately 50% with respect 
to the same calculation with pseudovector coupling. This is not, per se, ground enough to 
preclude the use of pseudoscalar coupling. However, for the many-body problem is more 
efficient [|l[] to work with a Lagrangian with pseudovector coupling and nocTvrvr term because 



in this form the soft-pion limits appear naturally in terms of just pions and nucleons (without 
the necessity of subtle cancellations with terms containing the a). 

One point to notice is the dependence of the MEC results on the value of g^. Calculations 
in finite nuclei within QHD (see, for example, Refs. and [^) usually employ = 65.23 
(the value used in this work) instead of g^ = 36.76 as calculated from the decay of the 
p meson. The first value is obtained fitting the symmetry energy of the empirical mass 
formula; this is a many-body argument which supports its use over that obtained from p 
decay. However, aside from this, there are no further arguments which are strong enough 
to discriminate them; this ambiguity translates into a downwards 40% uncertainty in the 
final MEC results. It is relevant to comment that the relativistic one-body calculation is 
not sensitive to the value of the p coupling constant. 



A final word about core polarization. In two previous works |T2| , p!3| we investigated the 
isovector weak response of the core to the external perturbation in the approximation of 
nuclear matter and in linear response theory summing the ring diagrams to all orders in 



the random phase approximation (RPA) — the so-called backflow correction ||25|j26|| . We 
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found that the effects of this core polarization (of positive and negative energy nucleons) 
were sizable for large values of the momentum transfer (for example those typical of the 
muon capture processes, (f ~ 0.9m^) but could be neglected for the small momentum 
transfers involved in nuclear (5 decay. Thus, we do not anticipate any dramatic increase in 
the calculated 6a beyond those presented in this work. 

Summarizing, we have studied the contribution of two-body (MEC) processes to axial 
charge transitions for nuclei in the lead, tin and oxygen regions. We conducted calculations 
in Dirac-Hartree approximation (the Walecka model) and in the relativistic Hartree approx- 
imation (where the full one-nucleon-loop effects are included.) Along the paper we discussed 
the sensitivities of the calculation to the change in some of the couplings and of the param- 
eters involved. With the choices described above the results we obtained indicate that one- 
and two-body processes enhance the matrix elements of the axial-charge operator by some 
(100±20) % in all three regions studied. This result agrees with the enhancement of the 
axial-coupling constant needed to understand theoretically a fit to eighteen first-forbidden 
13 decay transitions in the lead region. From the results of this calculation the presence of a 
heavy scalar meson to explain a physical effect gets further support. 
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FIGURES 

FIG. 1. a) and b): Pair graphs representing the exchange of a scalar meson. Diagrams c) and 
d) are those evaluated in relativistic mean- field theories (one body). 

FIG. 2. A "time-ordered" diagram involving an intermediate nucleon-antinucleon pair gener- 
ated by the presence of the scalar field. This diagram is included in the relativistic meson-exchange 
calculation. 

FIG. 3. Tree-level diagrams to be considered in the calculation of two-body processes: a) pvr, 
b) pTTTT, c) An, and d) Atttt. 

FIG. 4. a) ttNN, b) pNN and c) pmr vertices derived from the Lagrangian of QHD-II 
[Eq. in the text]. 

FIG. 5. The weak tts and irpe vertices needed for the calculation and not present explicitly in 
the QHD-II Lagrangian. 
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TABLES 



TABLE L Parameters used for the calculations presented in the text. 





n2 
ycrN 


rria 




9l 


3l 


M*/M 


MFT 


109.6 


520 


190.4 


178.4 


65.23 


0.54 


RHA 


54.3 


458 


102.8 


178.4 


65.23 


0.73 


TABLE IL 


Results for the 


: 2^9/2 


TT : I/19/2 


transition in ^'^^Pb - 


^ 209g^^ First row: 


results 



for the axial-charge operator. All the impulse approximation (lA) results, non-relativistic (NR), 
Dirac-Hartree (DH) and relativistic Hartree (RHA) are from reference [|l2|. The contribution of 
pion-exchange currents is summed in the columns labeled I A + MEC . The total increment 5a 
is given for the two relativistic calculations. Second and third rows: same as above but for the 
longitudinal {L\) and transverse electric {TY '^^) operators. 





NR 


DH 


RHA 


rpVel 


lA 
1.12x10^1 
5.61x10-3 
7.93x10-3 


lA lA+MEC Sa 
1.98x10-1 3.12x10-1 1.77 
5.77x10-3 
8.16x10-3 


lA lA+MEC 5a 
1.58x10-1 2.48x10-1 1.2 
5.69x10-3 
8.05x10-3 


TABLE III. Same as Table for the transition 207x1(1/2+) ^ 209pi3(i/2-) 




NR 


DH 


RHA 


^0 

LX 

fVel 


lA 

-0.75x10-1 
2.01x10-2 
2.85x10-2 


lA lA+MEC 5a 
-1.15x10-1 -1.60x10-1 1.12 
3.33x10-2 
4.71x10-2 


lA lA+MEC 5a 
-0.97x10-1 -1.37x10-1 0.81 
2.71x10-2 
3.83x10-2 
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TABLE IV. Same as Table for the transition ^33Sn(7/2-) ^ ^33313(7/2+) 





NR 


DH 


RHA 


fVel 


lA 
0.98x10-^ 
0.66x10-2 
0.93x10-2 


lA 
1.69x10-1 
0.81x10-2 
1.15x10-2 


lA+MEC 
2.22x10-1 


Sa 
1.26 


lA 
1.42x10-1 
0.78x10-2 
1.10x10-2 


lA+MEC 
1.70x10-1 


Sa 
0.73 



TABLE V. Same as Table H for the transition ^'5N(0") ^'^0(0+) except that the one-body 
longitudinal (LY) and transverse electric (Tj*^^') contributions have been substituted by the J = 3/2 
contribution from the transition due to the second term in Eq. ( |3.6| ) 





NR 


DH 


RHA 


C^j=3/2 


lA 
3.95x10-2 
-1.31x10-1 


lA lA+MEC Sa 
5.97x10-2 7.90x10-2 1.35 
-1.64x10-1 -1.78x10-1 0.36 


lA lA+MEC Sa 
5.86x10-2 7.51x10-2 0.90 
-1.56x10-1 -1.68x10-1 0.28 



TABLE VL Comparison of the MEC results using pseudovector (PV) coupling and pseu- 
doscalar plus the amr term for the axial charge operator. The results correspond to the 
u : 2(79/2 vr : I/19/2 transition in ^'^^Pb — > ^'^^Bi and the labels are the same as in Table 
The second row shows the enhancement with respect to the NR calculation. 



NR 


DH 


RHA 


lA 


lA 


MEC{PV) 


MEC{PS) 


lA 


MEC{PV) 


MEC{PS) 


-1.23x10"^ 


-1.98x10^1 


-1.14x10"^ 


-2.63x10"^ 


-1.57x10"^ 


-8.97x10-2 


-2.14x10"^ 




0.60 


0.90 


2.14 


0.28 


0.72 


1.72 
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